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We study the magnetic properties of Gadolinium- Yttrium Iron Garnet (Gd2,Y3_a;Fe50i2, x = 
3, 1.8) ferrite ceramics. The complex initial permeability is measured in the temperature range 2 K 
to 295 K at frequency of 1 kHz, and in the frequency range 100 Hz to 200 MHz at temperatures 4 K, 
77 K, and 295 K. The magnetic viscosity-induced imaginary part of the permeability is observed at 
low frequencies. Measurements of the magnetization noise are made at 4 K. Using the fluctuation- 
dissipation theorem, we find that the observed magnetization fluctuations are consistent with our 
measurements of the low-frequency imaginary part of permeability. We discuss some implications 
for proposed precision measurements as well as other possible applications. 



I. INTRODUCTION 

Rare-earth garnets have generated sustained interest 
in various soft-ferrite applications, and, more recently, in 

' more fundamental research. These materials combine at- 
tractive magnetic properties (large permeability, low loss 
angle) with high resistivity, which leads to their exten- 
sive use in microwave applications [l|, Q • Most of their 
applications up to now have been at room temperature, 
although there are some more recent suggestions for low- 
temperature applications also 0]. Gadolinium iron gar- 

' net (GdIG), for example, has found use in adiabatic de- 
magnetization refrigeration (ADR) in satellites [3| ■ From 
the perspective of more fundamental measurements, rare- 
earth garnets have generated great interest as attrac- 
tive materials for experiments searching for the parity 
and time-reversal invariance- violating permanent electric 
dipole moment of the electron [E S S B 9 > as well as 
for measurements of nuclear anapole moments llOll , and 
possibly measurements of the Weinberg angle [11| . 

The first evidence that GdIG retains its large initial 
permeability when it is cooled to cryogenic temperatures 
came from the study by Pascard [I4I of the real part of 
the magnetic susceptibility versus temperature, spanning 

[ the range from 4 K to 564 K (the Curie temperature of 
GdIG). The complex permeability ^ = /i' — ifi" of GdIG 

I was measured earlier at radio-frequencies at three tem- 
peratures: 300 K, 195 K, and 77 K [l^. Other properties 
have also been studied at low temperature, including the 
spontaneous and saturation magnetizations [l4l.fl5l , crys- 

' tal anisotropy [l(|, and ferrimagnetic order [l7l|7 How- 
ever, to our knowledge, there have been no studies of 
the full complex initial permeability of Gd-containing fer- 
rites in the temperature range between 4 K and 300 K, 
and at frequencies up to 1 GHz. Such a study is neces- 
sary to evaluate the suitability of these materials for low- 
temperature applications and precision measurements. 

An important emerging application of soft fcrrites 
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is magnetic shielding. Sensitive magnetic field mea- 
surements have to be performed inside a set of per- 
meable shields, and are often limited by the mag- 
netic field noise generated by Johnson currents in the 
shields themselves [3. Constructing the shields out of 
high-resistivity soft ferrite materials greatly reduces this 
source of magnetic noise at frequencies above approxi- 
mately 50 Hz 19]. At lower frequencies, however, finite 
imaginary part of the permeability generates extra mag- 
netic noise with 1/f power spectrurn, i n ag reement with 
the fluctuation-dissipation theorem [l^, [23| . 

This magnetic noise is related to magnetic viscosity, 
which arises when the magnetization of a sample is de- 
layed following application of a magnetic field [IJ]. Al- 
though initially seen in ferrous materials, this effect has 
recently been found in nanowires and other nanostuc- 
tures [23, ■ Anomalous magnetic viscosity effects have 
also been seen in various ferrite materials [24| . Magnetic 
viscosity is typically measured by looking at the repsonse 
of the magnetization after a sudden change in the applied 
field; however, it is also characterized by a small /i" that 
is independent of frequency for a given magnetic field 
strength and temperature [2^ [23] ■ As shown in the 
above references, the temperature and field dependences 
of this effect can be quite complex. 

In this work we study the magnetic properties of mixed 
Gadolinium and Yttrium iron garnets with the chemical 
formula Gda;Y3_a;Fe50i2. Here, x quantifies the ratio of 
Gadolinium to Ytrrium present within the garnet. We 
present measurements of both real and imaginary parts 
of the complex initial susceptibility in the temperature 
range of 4 K to 300 K, and in the frequency range of 
100 Hz to 200 MHz. We also measure the magnetization 
noise of these materials, and confirm the consistency of 
our results with the fiuctuation-dissipation theorem. 



II. EXPERIMENTAL SETUP 

The inductance of a tightly wound coil on a toroidal 
core is proportional to the permeability fi of the core 
material. The presence of a complex permeabililty fi = 
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FIG. 1: Equivalent circuit for an inductor with a dissipa- 
tive core (see text) and circuit schematic for the SQUID flux 
pickup loop, described in Section HVl 



fi' — ifi" creates a complex inductance L = L' — iL" . The 
impedance for such an toroidal inductor is Zt = iujL' + 
LuL", which is equivalent to a resistance of i? = ujL" 
in series with an inductor L', as shown in Fig. [TJ We 
deduce the real and imaginary parts of the permeability 
from measurements of the complex impedance of such a 
toroid. 

Two toroidal samples: one pure GdIG (Gd3Fe50i2), 
and one mixed GdYIG (Gdi.8Yi.2Fe50i2), were pur- 
chased from Pacific Ceramics. These were manufactured 
using their standard process. Powder was granulated and 
screened, using 99.99% pure Gd, 99.999% pure Y, and 
99.6% pure Fe oxides. The typical grain size was 5-10 
microns. The samples were pressed and fired at approxi- 
mately 1400°C for 15 hours. The toroids were made with 
square cross sections, with inner radius of 0.5 cm, outer 
radius of 1.5 cm, and a height of 1 cm. 

To measure the real part of the permeability, /i', at 
frequencies less than 1 MHz, the toroidal samples were 
wrapped with approximately 30 turns of wire, and the 
resulting self-inductance was measured with a QuadTech 
Series 1200 LCR meter. This method does not give reli- 
able measurements for the imaginary part of the perme- 
ability, however, since the ohmic resistance of the wind- 
ing dominates the loss due to fi". For measurements of 
n" at frequencies less than 1 MHz, the toroidal samples 
were wrapped with a primary winding of approximately 
35 turns, and a secondary winding of approximately 80 
turns. Litz wire was used to reduce skin effects. The 
complex mutual inductance of the two windings was then 
measured by a four-point measurement with the LCR 
meter. Since the input impedance of the voltage inputs 
on the meter is large (> 10 Mfl), the current flowing 
in the secondary winding is very small, thus the ohmic 
resistance of the secondary does not affect the measure- 
ment of the imaginary part of the susceptibility (on the 
level of a fraction of one percent). The inter- winding ca- 
pacitance, however, limits this technique to frequencies 
below 100 kHz, at higher frequencies the small impedance 
of the shunt capacitor results in a substantial current in 
the secondary winding. 

A Hewlett-Packard series 4100A impedance analyzer 
was used for measurements in the range of 1 MHz to 



200 MHz. The sample was placed in a toroidal-shaped 
copper housing to reduce the number of turns around 
the toroid to one. This minimized the electrical length of 
the wrapping such that it was always less than the wave- 
length. However, it introduced capacitive effects which 
prevented reliable measurements at frequencies above 
200 MHz. The self inductance of the single turn induc- 
tor was measured with the impedance analyzer. Ohmic 
resistance effects lead to errors on the order of 1%, since 
the series impedance due to /i" dominates at such high 
frequencies. 

For all the inductance measurements, the strength 
of the applied magnetic field was no more than Hq 
3 A/m, estimated from the current flowing in the toroid 
windings. This is much smaller than the ferrite's satu- 
ration field {Hq ^ 10^ A/m) pj]], which ensures that we 
always measure the initial susceptibility of the material. 

Low-temperature measurements were performed with 
the samples mounted in a G-10 holder inside a Janis 
model lOCNDT cryostat. Temperature was measured 
with a LakeShore model DT-670C-SD silicon diode tem- 
perature sensor, read out by a LakeShore Model 325 
temperature controller, with accuracy of 0.1 K. At 4 K, 
temperature-dependent effects on the calibration of the 
QuadTech LCR meter were corrected by making a mea- 
surement of tan (5 = fJ." / n' using a SQUID magnetome- 
ter. Multiple-turn coils were wrapped around toroidal 
non-magnetic (GIO) and GdIG samples and the coils 
were connected in series. A single-turn superconduct- 
ing pickup loop was wrapped about each sample and 
each pickup loop was connected to a SQUID. An sinu- 
soidal current was passed through the coils, and the phase 
shift between the non-mangetic and GdIG samples was 
recorded. The tan 5 recorded by the SQUIDs was then 
used to correct the phase-shift error in the calibration of 
the QuadTech LCR meter at 4 K. More details on the 
SQUID setup can be found in Sec. IV] and Ref. For 
temperature dependent measurments, temperature vari- 
ation was accomplished by removing the cryogens and 
letting the cryostat warm up, while a control computer 
read out the temperature controller and the LCR me- 
ter. Temperature dependent measurements were typi- 
cally separated by 0.1 K, and the total warming time 
was approximately 12 hours. 



III. INITIAL SUSCEPTIBILITY 

The measured real part of permeability, /i', as a func- 
tion of temperature is shown in Fig. [2] for GdIG and 
GdYIG ferrites. It is evident that both these ferrites 
retain their large permeability down to 4 K. The mag- 
netic properties of mixed Gadolinium- Yttrium iron gar- 
nets are well described by a model with three magnetic 
sublattices: the c-sites containing Gd'^"^ ions, the a-sites 
containing Fe'^+ ions, and the d-sites containing Fe'^^ 
ions [29]. The Y'^+ ions are not magnetic, they substitute 
for Gd'^^. At the compensation temperature, the magne- 
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FIG. 2: The real part of the initial permeability of GdIG 
and GdYIG versus temperature. The Globus model (Eq. [1]) 
described in the text is shown in fine dashed red as the best fit 
to the data between 65 K and 150 K. All data were obtained 
at a frequency of 1 kHz and H « 3 A/m applied. 



tizations of the three sublattices compensate each other, 
so that the net magnetization vanishes, and the perme- 
ability approaches unity js^]. As evident from Fig. [21 
the compensation temperature for Gdi.8Yi.2Fe50i2 is 
155 K, in good agreement with measurements in Ref. [l7[. 
The c omp ensation temperature for pure GdIG is 290 
K [ll,|IJ[3l|, which is not visible in the plot. 



As claimed in Ref. Il^ 
bulging due to Globus [32 



the model of domain wall 
33[ provides a satisfactory ex- 
planation of the temperature dependence of /x' for GdIG 
{x — 3.0). In this model 
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where Ms is the saturation magnetization, Ki is the mag- 
netocrystalline anisotropy. Dm is the mean grain size, 
and / is a constant, with dimensions of length, that de- 
pends on the properties of the material. Taking the data 
for Ms and Ki from Refs. [3, [l^ and fitting in the 
least squares sense the constant I yields the theoretical 
fit in Fig. [21 There is good agreement with our data be- 
tween 65 K and 200 K. Lack of agreement below 50 K, 
however, su gges ts that the single-crystal samples used 
in Refs. [13,13] do not have the same low-temperature 
properties as the pressed ceramic used for the present 
work. 

Our measurements indicate the maximum permeabil- 
ity of 112 at 47 K for GdIG. Susceptibility data from 
Ref. [l2], however, shows a peak in the GdIG susceptibil- 
ity of ii' ^ 70 at 56 K. This discrepancy can be accounted 
for within the model ([T]) by the difference in grain size. 
Dm, of the ceramic samples: the samples used in Ref. (T2j 
have Dm 4 ^m, whereas our samples have Dm ~ 7 /im. 



FIG. 3: The real and imaginary parts of permeability of GdIG 
versus applied field H, averaged over the volume of the sam- 
ple. This data was obtained at a frequency of 1 kHz and 
temperature of 4 K. 



At fixed temperature and frequency, the complex sus- 
ceptibility depends on the applied field H. The mea- 
surements of fj,' and fj," for various applied fields H are 
shown in Fig. [3l We always observed linear dependence 
of the magnetic susceptibility on H, implying that the 
magnetization is given by M = xH + ^aH'^ . This non- 
linear dependence of magnetization on applied field leads 
to the low-field Rayleigh hysteresis loop [3^]. Our four- 
point measurements of ^" at low frequency (/ < 100 
kHz) were always made at several applied fields H and 
linearly extrapolated to zero, as shown in Fig. [31 

The measurements of the complex susceptibility as a 
function of frequency at 4 K, 77 K, and 295 K are shown 
in Fig. m No data are shown for pure GdIG at 295 K, 
because this is close to its compensation temperature, so 
the sample was very nearly non-magnetic at this temper- 
ature. It is evident from the data that near 20 MHz the 
permeability of both samples exhibits a resonance. This 
"natural resonance" occurs at the frequency = 'jH^, 
where 7 is the gyromagnetic ratio, and cx Ki/Mg is 
the crystalline anisotropy field [U, [sH] ■ Using Eq. ([T]) we 
can re- write the resonance frequency as: 



Ms 



1 



(2) 



Using the data for saturation magnetization in Ref. [lj| , 
and our data for /i' — 1, it is apparent that as the tem- 
perature decreases, Ms grows faster than — 1, thus the 
natural resonance should shift to higher frequency with 
lower temperature, which is indeed the case as can be 
seen from the data in Fig. [H A similar trend is observed 
in Ref. pi]. 

In spite of the shift of the natural resonance to higher 
frequency, the low-frequency imaginary part of perme- 
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FIG. 4: The real and imaginary parts of permeability of GdIG 
and GdYIG versus frequency. Since /i' ~ 1 for GdIG at room 
temperature, no data is shown. All data were obtained with 
an applied field of // ~ 3 A/m. Low frequency (/ < 100 
kHz) ^" data were obtained by a four point measurement and 
extrapolated to zero applied field, as described in the text. 



ability increases as the temperature is lowered. This, as 
well as inspection of the data in Fig. 21 shows that the 
frequency dependence of the permeabihty is not as sim- 
ple as the models suggested in Ref. ^6], for example. It 
appears that at low frequencies (/ < 100 kHz) the mag- 
nitude of /i" is dominated not by the wing of the natural 
resonance, but by some frequency-independent dissipa- 
tion mechanism. This is known as magnetic viscosity. 

The presence of a non-zero fi" is indicative of a phase 
shift of magnetization with respect to the applied field H. 
The phase lag, or loss angle, S is given by tan 6 = fJ." / n' ■ 
For a given temperature, and at frequencies much lower 
than the natural resonance, the loss angle is empirically 
parameterized by the lag equation (st} : 



„ tan (5 

27r = c + aH + ef, 



(3) 



where c, a, and e are parameters, and / is the frequency 
of the applied field H, which is much less than the satura- 
tion field. The last term accounts for eddy current effects. 
However, in the present case, eddy current losses are neg- 
ligible since our ferrite samples are very good insulators. 



TABLE I: Parameters of the loss equation (|3]) for GdIG and 
mixed GdYIG at various temperatures. Values were compiled 
by least squares fit to all available data less than 10 kHz. 
Error bars are extracted from the fit (1(j). 
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The second term describes hysteresis effects described 
above, thus the parameter a is a measure of non-linearity 
in the system. The parameter c is a constant allowing for 
frequency-independent loss at zero applied field, this is 
the magnetic viscosity. Using the low-frequency measure- 
ments of /i" versus frequency and applied field H , values 
of a and c were extracted at 4 K, 77 K, and 295 K. These 
values are shown in Table HI Once again, we do not give 
any results for pure GdIG {x = 3) at 295 K, since this is 
very close to its compensation temperature. 

Magnetic viscosity is characterized by time dependent 
change in the magnetization after a change in the ap- 
plied magnetic field. When a small magnetic field is ap- 
plied to a ferromagnet, magnetization changes primarily 
due to domain wall movement. As it moves, the domain 
wall traverses a complex potential energy landscape [s^ . 
The time scale of the thermally-driven movement across 
a potential energy barrier of height Ed is given by an 
Arrhenius-type equation: t = TQex.'p{Ed/kBT) [35'|. In 
a macroscopic sample, there are multiple domain walls, 
each seeing some distribution of such activation energies. 
If this distribution is flat over a wide range of potential 
barrier heights, then /x" is independent of frequency, and 
the magnetic viscosity term c a ppe ars in the expression 
dSl) for the lag angle [llllElUlsil. 

The temperature dependence of magnetic viscosity of 
our GdYIG samples can be inferred from Fig.[5l where ^" 
is plotted as a function of temperature. There are several 
competing effects that lead to a complicated temperature 
dependence [25|. On the one hand, at low temperature, 
the domain wall relaxation times grow longer, increasing 
the magnitude of /i". On the other hand, as the temper- 
ature is lowered even further, some of the domain walls 
become pinned, no longer contributing to the initial sus- 
ceptibility, decreasing the magnitudes of both /i' and /i", 
as seen in the x = 3.0 curve in Figs. and [21 



IV. MAGNETIZATION NOISE 

The magnitude of the magnetization noise is the key 
property of the GdYIG garnet ferrites that determines 
the feasibility of their application in magnetic shielding as 
well as in experiments searching for parity and/or time- 
reversal symmetry violations. We measured the magne- 
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FIG. 5: Measurement of /i" versus temperature for GdIG and 
GdYIG. All data were obtained at a frequency of 1 kHz and 
// ~ 3 A/m applied. Data below 25 K was unreliable due to 
significant temperature dependent effects on the QuadTech 
LCR meter's calibration. 
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FIG. 6: Magnetization noise of GdIG at 4.2 K. The blue 
line is the result of using the fluctuation-dissipation theorem, 
Eq. ([5]), with the value for the imaginary part of permeability 
fi" — 1.1, taken from measurements shown in Fig |4l The 
SQUID white noise level corresponds to 10 fT/v'Hz on this 
plot. 



tization noise of our samples and compared the results 
with our data for the complex initial permeability using 
the fluctuation-dissipation theorem. 

The noise measurements were made at 4 K inside the 
cryogenic dewar described in Section [ill A Quantum De- 
sign model 50 superconducting quantum interference de- 
vice (SQUID) magnetometer was connected to the sam- 
ple under study by a niobium wire pickup loop. The one- 
turn loop was wound around the square cross-section of 
the toroidal sample, so that the total flux pickup area 
was 1 cm^. The SQUID control and readout was per- 
formed by the control unit (model 5000), connected to 
the SQUID via a 4-meter MicroPREAMP cable. The 
SQUID and the sample were mounted on G-10 holders, 
and enclosed in superconducting shielding, made of 0.1- 
mm thick Pb foil, glued to the inner surfaces of two G- 
10 cylinders. The magnetic field shielding factor of our 
setup is greater than 10®. Some more details on our noise 
measurement setup are given in Ref. [28j . 

The measured magnetization noise amplitude spec- 
trum for a pure GdIG sample is shown in Fig. [6l The 
noise data for the mixed GdYIG sample is within 40% of 
the pure GdIG noise; we do not show it, since the two 
curves would lie on top of each other, given the plot's 
logarithmic scale. The measurement bandwidth of 1 kHz 
is set by the output analog filter to prevent aliasing with 
our sampling rate of 12 kHz. The SQUID magnetometer 
intrinsic sensitivity in this frequency range is 7 fT/VHz; 
within our bandwidth, the magnetization noise is always 
above the magnetometer sensitivity. 

According to the fluctuation-dissipation theorem, any 
dissipative material exhibits thermal noise. We measured 
the dissipation of our ferrite samples in Section [Till this 
is given by the imaginary part /i" of the complex perme- 



ability. We can therefore check that our measurements 
of the magnetization noise and the dissipation are consis- 
tent using the fluctuation-dissipation theorem. Consider 
the circuit shown in Fig.[ll where the complex impedance 
Zt of the pickup loop around the permeable sample is in 
series with the input inductance Li„ of the SQUID. The 
real part of Zt, ujL", is a source of voltage noise, whose 
spectral density is given by the Nyquist theorem: 

(SV)^ = ^AksTLoL", (4) 

where ks is Boltzmann's constant, T is temperature, 
Lu — 27r/ is angular frequency, and L" is the imaginary 
part of the pickup loop inductance. The current noise 
in the circuit can be obtained by dividing Eq. (j4|) by the 
total impedance, Zt + a;Lj„ w uj{L' + L^n). This current 
couples to the SQUID through the mutual inductance 
Msg. The resulting flux through the SQUID is the same 
as would be created by a lossless pickup loop around the 
sample with magnetization noise given by: 

fioA y 27r/ 

where A = 1.0 cm^ is the area of the pickup loop and 
fiQ is the permeability of free space. The dependence of 
the single turn pickup loop's inductance on x = ^ 1 
was found empirically and is described well by a linear 
expansion, e.g. L = Lq + x^i, where Lq and Li are 
the expansion coeffecients. Since these cofFecients are 
real, the imaginary part of the inductance is given by 

Using Eq. ^ and our empirical determination of L" , 
we can compare our data for /i" of GdIG at 4 K with 
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our measurement of its magnetization noise. From the 
GdIG data shown in Fig. 0] we can deduce, with good ac- 
curacy, a constant /i" = 1.1 for frequencies below 1 kHz. 
With this value, Eq. (O gives a prediction for the mag- 
netization noise, shown as the blue line in Fig. [G] Note 
the l/\ff dependence characteristic of magnetic viscos- 
ity. The agreement with experimental noise measure- 
ments is remarkable, given that fx" was measured with 
an entirely different method. 

V. CONCLUSIONS 

We have measured the full complex permeability of two 
GdYIG ceramics in the temperature range 2 K to 295 K 
and confirmed that the real part of the permeability re- 
mains large (78 and 52 respectively) at 4 K. The Globus 
magnetization model fits the results well between 65 K 
and 200 K. We also measured the full complex perme- 
ability as a function of frequency in the range 100 Hz to 
200 MHz at temperatures of 4 K, 77 K, and 295 K. We 
observed the natural resonances near 20 MHz, and de- 
tected the frequency-independent imaginary permeabil- 
ity due to magnetic viscosity effects for frequencies less 
than 100 kHz: /x" = 1.1 for GdlG and /z" = 0.8 for 
GdYIG at 4 K. To our knowledge, such effects have not 
been studied for rare-earth ferrites in such a wide temper- 
ature range. Using the fluctuation-dissipation theorem, 
these results were compared with direct measurements of 



the magnetization noise, with very good agreement. 

Our noise measurements at 4 K show that the intrinsic 
1 / -/J magnetization noise in GdIG and GdYIG is larger 
than the SQUID magnetometer noise for frequencies be- 
low approximately 100 kHz. This implies that precision 
measurements involving these materials either have to be 
carried out at higher frequencies, to avoid this noise, or 
have to make use of the geometry-dependent demagne- 
tizing fields to reduce the coupling of the noise to the 
measurement device [28^. 

A possible application of the materials studied in the 
present work is to magnetic shielding of high-Tc SQUIDs 
at 77 K. For a cylindrical shield of radius 5 cm and thick- 
ness 1 cm (the dimensions in Ref. [l3)j made of GdIG, 
the shielding factor is 6. The magnetic field noise on the 
axis can be predicted from our measurement of ^" = 0.5 
at 77 K, resulting in 5B = 10 fr/\/Hz at 1 Hz. This is at 
the level of a high-Tc SQUID's sensitivity. Such a shield 
can be used inside a set of /i-metal shields to screen ex- 
ternal magnetic fields and the Johnson noise due to the 
innermost layer of the /x-metal shields, without generat- 
ing excessive Johnson noise itself. 
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